Cellular heterogeneity and dynamics of the human uterus in healthy premenopausal women.

The human uterus is a complex and dynamic organ whose lining grows, remodels, and regenerates in every menstrual cycle or upon tissue damage. Here we applied single-cell RNA sequencing to profile more the 50,000 uterine cells from both the endometrium and myometrium of 5 healthy premenopausal individuals, and jointly analyzed the data with a previously published dataset from 15 subjects. The resulting normal uterus cell atlas contains more than 167K cells representing the lymphatic endothelium, blood endothelium, stromal, ciliated epithelium, unciliated epithelium, and immune cell populations. Focused analyses within each major cell type and comparisons with subtype labels from prior studies allowed us to document supporting evidence, resolve naming conflicts, and to propose a consensus annotation system of 39 subtypes. We release their gene expression centroids, differentially expressed genes, and mRNA patterns of literature-based markers as a shared community resource. We find many subtypes show dynamic changes over different phases of the cycle and identify multiple potential progenitor cells: compartment-wide progenitors for each major cell type, transitional cells that are upstream of other subtypes, and potential cross-lineage multipotent stromal progenitors that may be capable of replenishing the epithelial, stromal, and endothelial compartments. When compared to the healthy premenopausal samples, a postpartum and a postmenopausal uterus sample revealed substantially altered tissue composition, involving the rise or fall of stromal, endothelial, and immune cells. The cell taxonomy and molecular markers we report here are expected to inform studies of both basic biology of uterine function and its disorders. SIGNIFICANCE: We present single-cell RNA sequencing data from seven individuals (five healthy pre-menopausal women, one post-menopausal woman, and one postpartum) and perform an integrated analysis of this data alongside 15 previously published scRNA-seq datasets. We identified 39 distinct cell subtypes across four major cell types in the uterus. By using RNA velocity analysis and centroid-centroid comparisons we identify multiple computationally predicted progenitor populations for each of the major cell compartments, as well as potential cross-compartment, multi-potent progenitors. While the function and interactions of these cell populations remain to be validated through future experiments, the markers and their "dual characteristics" that we describe will serve as a rich resource to the scientific community. Importantly, we address a significant challenge in the field: reconciling multiple uterine cell taxonomies being proposed. To achieve this, we focused on integrating historical and contemporary knowledge across multiple studies. By providing detailed evidence used for cell classification we lay the groundwork for establishing a stable, consensus cell atlas of the human uterus.

Cellular atlas of the human ovary using morphologically guided spatial transcriptomics and single-cell sequencing.

The reproductive and endocrine functions of the ovary involve spatially defined interactions among specialized cell populations. Despite the ovary's importance in fertility and endocrine health, functional attributes of ovarian cells are largely uncharacterized. Here, we profiled >18,000 genes in 257 regions from the ovaries of two premenopausal donors to examine the functional units in the ovary. We also generated single-cell RNA sequencing data for 21,198 cells from three additional donors and identified four major cell types and four immune cell subtypes. Custom selection of sampling areas revealed distinct gene activities for oocytes, theca, and granulosa cells. These data contributed panels of oocyte-, theca-, and granulosa-specific genes, thus expanding the knowledge of molecular programs driving follicle development. Serial samples around oocytes and across the cortex and medulla uncovered previously unappreciated variation of hormone and extracellular matrix remodeling activities. This combined spatial and single-cell atlas serves as a resource for future studies of rare cells and pathological states in the ovary.

A revamped rat reference genome improves the discovery of genetic diversity in laboratory rats.

The seventh iteration of the reference genome assembly for Rattus norvegicus-mRatBN7.2-corrects numerous misplaced segments and reduces base-level errors by approximately 9-fold and increases contiguity by 290-fold compared with its predecessor. Gene annotations are now more complete, improving the mapping precision of genomic, transcriptomic, and proteomics datasets. We jointly analyzed 163 short-read whole-genome sequencing datasets representing 120 laboratory rat strains and substrains using mRatBN7.2. We defined ∼20.0 million sequence variations, of which 18,700 are predicted to potentially impact the function of 6,677 genes. We also generated a new rat genetic map from 1,893 heterogeneous stock rats and annotated transcription start sites and alternative polyadenylation sites. The mRatBN7.2 assembly, along with the extensive analysis of genomic variations among rat strains, enhances our understanding of the rat genome, providing researchers with an expanded resource for studies involving rats.

Deriving early single-rosette brain organoids from human pluripotent stem cells.

Brain organoid methods are complicated by multiple rosette structures and morphological variability. We have developed a human brain organoid technique that generates self-organizing, single-rosette cortical organoids (SOSR-COs) with reproducible size and structure at early timepoints. Rather than patterning a 3-dimensional embryoid body, we initiate brain organoid formation from a 2-dimensional monolayer of human pluripotent stem cells patterned with small molecules into neuroepithelium and differentiated to cells of the developing dorsal cerebral cortex. This approach recapitulates the 2D to 3D developmental transition from neural plate to neural tube. Most monolayer fragments form spheres with a single central lumen. Over time, the SOSR-COs develop appropriate progenitor and cortical laminar cell types as shown by immunocytochemistry and single-cell RNA sequencing. At early time points, this method demonstrates robust structural phenotypes after chemical teratogen exposure or when modeling a genetic neurodevelopmental disorder, and should prove useful for studies of human brain development and disease modeling.

Sperm chromatin structure and reproductive fitness are altered by substitution of a single amino acid in mouse protamine 1.

Conventional dogma presumes that protamine-mediated DNA compaction in sperm is achieved by electrostatic interactions between DNA and the arginine-rich core of protamines. Phylogenetic analysis reveals several non-arginine residues conserved within, but not across species. The significance of these residues and their post-translational modifications are poorly understood. Here, we investigated the role of K49, a rodent-specific lysine residue in protamine 1 (P1) that is acetylated early in spermiogenesis and retained in sperm. In sperm, alanine substitution (P1(K49A)) decreases sperm motility and male fertility-defects that are not rescued by arginine substitution (P1(K49R)). In zygotes, P1(K49A) leads to premature male pronuclear decompaction, altered DNA replication, and embryonic arrest. In vitro, P1(K49A) decreases protamine-DNA binding and alters DNA compaction and decompaction kinetics. Hence, a single amino acid substitution outside the P1 arginine core is sufficient to profoundly alter protein function and developmental outcomes, suggesting that protamine non-arginine residues are essential for reproductive fitness.

A revamped rat reference genome improves the discovery of genetic diversity in laboratory rats.

The seventh iteration of the reference genome assembly for Rattus norvegicus-mRatBN7.2-corrects numerous misplaced segments and reduces base-level errors by approximately 9-fold and increases contiguity by 290-fold compared to its predecessor. Gene annotations are now more complete, significantly improving the mapping precision of genomic, transcriptomic, and proteomics data sets. We jointly analyzed 163 short-read whole genome sequencing datasets representing 120 laboratory rat strains and substrains using mRatBN7.2. We defined ~20.0 million sequence variations, of which 18.7 thousand are predicted to potentially impact the function of 6,677 genes. We also generated a new rat genetic map from 1,893 heterogeneous stock rats and annotated transcription start sites and alternative polyadenylation sites. The mRatBN7.2 assembly, along with the extensive analysis of genomic variations among rat strains, enhances our understanding of the rat genome, providing researchers with an expanded resource for studies involving rats.

The Legacy of Infectious Disease Exposure on the Genomic Diversity of Indigenous Southern Mexicans.

To characterize host risk factors for infectious disease in Mesoamerican populations, we interrogated 857,481 SNPs assayed using the Affymetrix 6.0 genotyping array for signatures of natural selection in immune response genes. We applied three statistical tests to identify signatures of natural selection: locus-specific branch length (LSBL), the cross-population extended haplotype homozygosity (XP-EHH), and the integrated haplotype score (iHS). Each of the haplotype tests (XP-EHH and iHS) were paired with LSBL and significance was determined at the 1% level. For the paired analyses, we identified 95 statistically significant windows for XP-EHH/LSBL and 63 statistically significant windows for iHS/LSBL. Among our top immune response loci, we found evidence of recent directional selection associated with the major histocompatibility complex (MHC) and the peroxisome proliferator-activated receptor gamma (PPAR-γ) signaling pathway. These findings illustrate that Mesoamerican populations' immunity has been shaped by exposure to infectious disease. As targets of selection, these variants are likely to encode phenotypes that manifest themselves physiologically and therefore may contribute to population-level variation in immune response. Our results shed light on past selective events influencing the host response to modern diseases, both pathogenic infection as well as autoimmune disorders.

Presence of ovarian stromal aberrations after cessation of testosterone therapy in a transgender mouse model†.

Some transmasculine individuals may be interested in pausing gender-affirming testosterone therapy and carrying a pregnancy. The ovarian impact of taking and pausing testosterone is not completely understood. The objective of this study was to utilize a mouse model mimicking transmasculine testosterone therapy to characterize the ovarian dynamics following testosterone cessation. We injected postpubertal 9-10-week-old female C57BL/6N mice once weekly with 0.9 mg of testosterone enanthate or a vehicle control for 6 weeks. All testosterone-treated mice stopped cycling and demonstrated persistent diestrus within 1 week of starting testosterone, while control mice cycled regularly. After 6 weeks of testosterone therapy, one group of testosterone-treated mice and age-matched vehicle-treated diestrus controls were sacrificed. Another group of testosterone-treated mice were maintained after stopping testosterone therapy and were sacrificed in diestrus four cycles after the resumption of cyclicity along with age-matched vehicle-treated controls. Ovarian histological analysis revealed stromal changes with clusters of large round cells in the post testosterone group as compared to both age-matched controls and mice at 6 weeks on testosterone. These clusters exhibited periodic acid-Schiff staining, which has been previously reported in multinucleated macrophages in aging mouse ovaries. Notably, many of these cells also demonstrated positive staining for macrophage markers CD68 and CD11b. Ovarian ribonucleic acid-sequencing found upregulation of immune pathways post testosterone as compared to age-matched controls and ovaries at 6 weeks on testosterone. Although functional significance remains unknown, further attention to the ovarian stroma may be relevant for transmasculine people interested in pausing testosterone to carry a pregnancy.

Intrinsic cardiorespiratory fitness modulates clinical and molecular response to caloric restriction.

OBJECTIVE: Caloric restriction (CR) is one extrinsic intervention that can improve metabolic health, and it shares many phenotypical parallels with intrinsic high cardiorespiratory fitness (CRF), including reduced adiposity, increased cardiometabolic health, and increased longevity. CRF is a highly heritable trait in humans and has been established in a genetic rat model selectively bred for high (HCR) and low (LCR) CRF, in which the HCR live longer and have reduced body weight compared to LCR. This study addresses whether the inherited high CRF phenotype occurs through similar mechanisms by which CR promotes health and longevity. METHODS: We compared HCR and LCR male rats fed ad libitum (AL) or calorically restricted (CR) for multiple physiological, metabolic, and molecular traits, including running capacity at 2, 8, and 12 months; per-hour metabolic cage activity over daily cycles at 6 and 12 months; and plasma lipidomics, liver and muscle transcriptomics, and body composition after 12 months of treatment. RESULTS: LCR-CR developed a physiological profile that mirrors the high-CRF phenotype in HCR-AL, including reduced adiposity and increased insulin sensitivity. HCR show higher spontaneous activity than LCR. Temporal modeling of hourly energy expenditure (EE) dynamics during the day, adjusted for body weight and hourly activity levels, suggest that CR has an EE-suppressing effect, and high-CRF has an EE-enhancing effect. Pathway analysis of gene transcripts indicates that HCR and LCR both show a response to CR that is similar in the muscle and different in the liver. CONCLUSIONS: CR provides LCR a health-associated positive effect on physiological parameters that strongly resemble HCR. Analysis of whole-body EE and transcriptomics suggests that HCR and LCR show line-dependent responses to CR that may be accreditable to difference in genetic makeup. The results do not preclude the possibility that CRF and CR pathways may converge.

Fibromuscular Dysplasia and Abdominal Aortic Aneurysms Are Dimorphic Sex-Specific Diseases With Shared Complex Genetic Architecture.

BACKGROUND: The risk of arterial diseases may be elevated among family members of individuals having multifocal fibromuscular dysplasia (FMD). We sought to investigate the risk of arterial diseases in families of individuals with FMD. METHODS: Family histories for 73 probands with FMD were obtained, which included an analysis of 463 total first-degree relatives focusing on FMD and related arterial disorders. A polygenic risk score for FMD (PRSFMD) was constructed from prior genome-wide association findings of 584 FMD cases and 7139 controls and evaluated for association with an abdominal aortic aneurysm (AAA) in a cohort of 9693 AAA cases and 294 049 controls. A previously published PRSAAA was also assessed among the FMD cases and controls. RESULTS: Of all first degree relatives of probands, 9.3% were diagnosed with FMD, aneurysms, and dissections. Aneurysmal disease occurred in 60.5% of affected relatives and 5.6% of all relatives. Among 227 female first-degree relatives of probands, 4.8% (11) had FMD, representing a relative risk (RR)FMD of 1.5 ([95% CI, 0.75-2.8]; P=0.19) compared with the estimated population prevalence of 3.3%, though not of statistical significance. Of all fathers of FMD probands, 11% had AAAs resulting in a RRAAA of 2.3 ([95% CI, 1.12-4.6]; P=0.014) compared with population estimates. The PRSFMD was found to be associated with an AAA (odds ratio, 1.03 [95% CI, 1.01-1.05]; P=2.6×10-3), and the PRSAAA was found to be associated with FMD (odds ratio, 1.53 [95% CI, 1.2-1.9]; P=9.0×10-5) as well. CONCLUSIONS: FMD and AAAs seem to be sex-dimorphic manifestations of a heritable arterial disease with a partially shared complex genetic architecture. Excess risk of having an AAA according to a family history of FMD may justify screening in family members of individuals having FMD.

Burden of Rare Genetic Variants in Spontaneous Coronary Artery Dissection With High-risk Features.

Importance: The emerging genetic basis of spontaneous coronary artery dissection (SCAD) has been defined as both partially complex and monogenic in some patients, involving variants predominantly in genes known to underlie vascular connective tissue diseases (CTDs). The effect of these genetic influences has not been defined in high-risk SCAD phenotypes, and the identification of a high-risk subgroup of individuals may help to guide clinical genetic evaluations of SCAD. Objective: To identify and quantify the burden of rare genetic variation in individuals with SCAD with high-risk clinical features. Design, Setting, and Participants: Whole-exome sequencing (WES) was performed for subsequent case-control association analyses and individual variant annotation among individuals with high-risk SCAD. Genetic variants were annotated for pathogenicity by in-silico analysis of genes previously defined by sequencing for vascular CTDs and/or SCAD, as well as genes prioritized by genome-wide association study (GWAS) and colocalization of arterial expression quantitative trait loci. Unbiased genome-wide association analysis of the WES data was performed by comparing aggregated variants in individuals with SCAD to healthy matched controls or the Genome Aggregation Database (gnomAD). This study was conducted at a tertiary care center. Individuals in the Canadian SCAD Registry genetics study with a high-risk SCAD phenotype were selected and defined as peripartum SCAD, recurrent SCAD, or SCAD in an individual with family history of arteriopathy. Main Outcomes and Measures: Burden of genetic variants defined by DNA sequencing in individuals with high-risk SCAD. Results: This study included a total of 336 participants (mean [SD] age, 53.0 [9.5] years; 301 female participants [90%]). Variants in vascular CTD genes were identified in 17.0% of individuals (16 of 94) with high-risk SCAD and were enriched (OR, 2.6; 95% CI, 1.6-4.2; P = 7.8 × 10-4) as compared with gnomAD, with leading significant signals in COL3A1 (OR, 13.4; 95% CI, 4.9-36.2; P = 2.8 × 10-4) and Loeys-Dietz syndrome genes (OR, 7.9; 95% CI, 2.9-21.2; P = 2.0 × 10-3). Variants in GWAS-prioritized genes, observed in 6.4% of individuals (6 of 94) with high-risk SCAD, were also enriched (OR, 3.6; 95% CI, 1.6-8.2; P = 7.4 × 10-3). Variants annotated as likely pathogenic or pathogenic occurred in 4 individuals, in the COL3A1, TGFBR2, and ADAMTSL4 genes. Genome-wide aggregated variant testing identified novel associations with peripartum SCAD. Conclusions and Relevance: In this genetic study, approximately 1 in 5 individuals with a high-risk SCAD phenotype harbored a rare genetic variant in genes currently implicated for SCAD. Genetic screening in this subgroup of individuals presenting with SCAD may be considered.

Decoding the Spermatogenesis Program: New Insights from Transcriptomic Analyses.

Spermatogenesis is a complex differentiation process coordinated spatiotemporally across and along seminiferous tubules. Cellular heterogeneity has made it challenging to obtain stage-specific molecular profiles of germ and somatic cells using bulk transcriptomic analyses. This has limited our ability to understand regulation of spermatogenesis and to integrate knowledge from model organisms to humans. The recent advancement of single-cell RNA-sequencing (scRNA-seq) technologies provides insights into the cell type diversity and molecular signatures in the testis. Fine-grained cell atlases of the testis contain both known and novel cell types and define the functional states along the germ cell developmental trajectory in many species. These atlases provide a reference system for integrated interspecies comparisons to discover mechanistic parallels and to enable future studies. Despite recent advances, we currently lack high-resolution data to probe germ cell-somatic cell interactions in the tissue environment, but the use of highly multiplexed spatial analysis technologies has begun to resolve this problem. Taken together, recent single-cell studies provide an improvedunderstanding of gametogenesis to examine underlying causes of infertility and enable the development of new therapeutic interventions.

Cellular heterogeneity of human fallopian tubes in normal and hydrosalpinx disease states identified using scRNA-seq.

Fallopian tube (FT) homeostasis requires dynamic regulation of heterogeneous cell populations and is disrupted in infertility and ovarian cancer. Here, we applied single-cell RNA-seq to profile 59,738 FT cells from four healthy, pre-menopausal subjects. The resulting cell atlas contains 12 major cell types representing epithelial, stromal, and immune compartments. Re-clustering of epithelial cells identified four ciliated and six non-ciliated secretory epithelial subtypes, two of which represent potential progenitor pools: one leading to mature secretory cells and the other contributing to either ciliated cells or one of the stromal cell types. To understand how FT cell numbers and states change in a disease state, we analyzed 17,798 cells from two hydrosalpinx samples and observed shifts in epithelial and stromal populations and cell-type-specific changes in extracellular matrix and TGF-ß signaling; this underscores fibrosis pathophysiology. This resource is expected to facilitate future studies aimed at expanding understanding of fallopian tube homeostasis in normal development and disease.

Molecular genetic evaluation of pediatric renovascular hypertension due to renal artery stenosis and abdominal aortic coarctation in neurofibromatosis type 1.

The etiology of renal artery stenosis (RAS) and abdominal aortic coarctation (AAC) causing the midaortic syndrome (MAS), often resulting in renovascular hypertension (RVH), remains ill-defined. Neurofibromatosis type 1 (NF-1) is frequently observed in children with RVH. Consecutive pediatric patients (N = 102) presenting with RVH secondary to RAS with and without concurrent AAC were prospectively enrolled in a clinical data base, and blood, saliva and operative tissue, when available, were collected. Among the 102 children, 13 were having a concurrent clinical diagnosis of NF-1 (12.5%). Whole exome sequencing was performed for germline variant detection, and RNA-Seq analysis of NF1, MAPK pathway genes and MCP1 levels were undertaken in five NF-1 stenotic renal arteries, as well as control renal and mesenteric arteries from children with no known vasculopathy or NF-1. In 11 unrelated children with sequencing data, 11 NF1 genetic variants were identified, of which 10 had not been reported in gnomAD. Histologic analysis of NF-1 RAS specimens consistently revealed intimal thickening, disruption of the internal elastic lamina and medial thinning. Analysis of transcript expression in arterial lesions documented an approximately 5-fold reduction in NF1 expression, confirming heterozygosity, MAPK pathway activation and increased MCP1 expression. In summary, NF-1-related RVH in children is rare but often severe and progressive and, as such, important to recognize. It is associated with histologic and molecular features consistent with an aggressive adverse vascular remodeling process. Further research is necessary to define the mechanisms underlying these findings.

Genetic investigation of fibromuscular dysplasia identifies risk loci and shared genetics with common cardiovascular diseases.

Fibromuscular dysplasia (FMD) is an arteriopathy associated with hypertension, stroke and myocardial infarction, affecting mostly women. We report results from the first genome-wide association meta-analysis of six studies including 1556 FMD cases and 7100 controls. We find an estimate of SNP-based heritability compatible with FMD having a polygenic basis, and report four robustly associated loci (PHACTR1, LRP1, ATP2B1, and LIMA1). Transcriptome-wide association analysis in arteries identifies one additional locus (SLC24A3). We characterize open chromatin in arterial primary cells and find that FMD associated variants are located in arterial-specific regulatory elements. Target genes are broadly involved in mechanisms related to actin cytoskeleton and intracellular calcium homeostasis, central to vascular contraction. We find significant genetic overlap between FMD and more common cardiovascular diseases and traits including blood pressure, migraine, intracranial aneurysm, and coronary artery disease.

Allelic Variants in Established Hypopituitarism Genes Expand Our Knowledge of the Phenotypic Spectrum.

We report four allelic variants (three novel) in three genes previously established as causal for hypopituitarism or related disorders. A novel homozygous variant in the growth hormone gene, GH1 c.171delT (p.Phe 57Leufs*43), was found in a male patient with severe isolated growth hormone deficiency (IGHD) born to consanguineous parents. A hemizygous SOX3 allelic variant (p.Met304Ile) was found in a male patient with IGHD and hypoplastic anterior pituitary. YASARA, a tool to evaluate protein stability, suggests that p.Met304Ile destabilizes the SOX3 protein (ΔΔG = 2.49 kcal/mol). A rare, heterozygous missense variant in the TALE homeobox protein gene, TGIF1 (c.268C>T:p.Arg90Cys) was found in a patient with combined pituitary hormone deficiency (CPHD), diabetes insipidus, and syndromic features of holoprosencephaly (HPE). This variant was previously reported in a patient with severe holoprosencephaly and shown to affect TGIF1 function. A novel heterozygous TGIF1 variant (c.82T>C:p.Ser28Pro) was identified in a patient with CPHD, pituitary aplasia and ectopic posterior lobe. Both TGIF1 variants have an autosomal dominant pattern of inheritance with incomplete penetrance. In conclusion, we have found allelic variants in three genes in hypopituitarism patients. We discuss these variants and associated patient phenotypes in relation to previously reported variants in these genes, expanding our knowledge of the phenotypic spectrum in patient populations.

High-throughput splicing assays identify missense and silent splice-disruptive POU1F1 variants underlying pituitary hormone deficiency.

Pituitary hormone deficiency occurs in ∼1:4,000 live births. Approximately 3% of the cases are due to mutations in the alpha isoform of POU1F1, a pituitary-specific transcriptional activator. We found four separate heterozygous missense variants in unrelated individuals with hypopituitarism that were predicted to affect a minor isoform, POU1F1 beta, which can act as a transcriptional repressor. These variants retain repressor activity, but they shift splicing to favor the expression of the beta isoform, resulting in dominant-negative loss of function. Using a high-throughput splicing reporter assay, we tested 1,070 single-nucleotide variants in POU1F1. We identified 96 splice-disruptive variants, including 14 synonymous variants. In separate cohorts, we found two additional synonymous variants nominated by this screen that co-segregate with hypopituitarism. This study underlines the importance of evaluating the impact of variants on splicing and provides a catalog for interpretation of variants of unknown significance in POU1F1.

TCF21+ mesenchymal cells contribute to testis somatic cell development, homeostasis, and regeneration in mice.

Testicular development and function rely on interactions between somatic cells and the germline, but similar to other organs, regenerative capacity declines in aging and disease. Whether the adult testis maintains a reserve progenitor population remains uncertain. Here, we characterize a recently identified mouse testis interstitial population expressing the transcription factor Tcf21. We found that TCF21lin cells are bipotential somatic progenitors present in fetal testis and ovary, maintain adult testis homeostasis during aging, and act as potential reserve somatic progenitors following injury. In vitro, TCF21lin cells are multipotent mesenchymal progenitors which form multiple somatic lineages including Leydig and myoid cells. Additionally, TCF21+ cells resemble resident fibroblast populations reported in other organs having roles in tissue homeostasis, fibrosis, and regeneration. Our findings reveal that the testis, like other organs, maintains multipotent mesenchymal progenitors that can be potentially leveraged in development of future therapies for hypoandrogenism and/or infertility.

The phenotypic spectrum associated with OTX2 mutations in humans.

Objective: The transcription factor OTX2 is implicated in ocular, craniofacial, and pituitary development. Design: We aimed to establish the contribution of OTX2 mutations in congenital hypopituitarism patients with/without eye abnormalities, study functional consequences, and establish OTX2 expression in the human brain, with a view to investigate the mechanism of action. Methods: We screened patients from the UK (n = 103), international centres (n = 24), and Brazil (n = 282); 145 were within the septo-optic dysplasia spectrum, and 264 had no eye phenotype. Transactivation ability of OTX2 variants was analysed in murine hypothalamic GT1-7 neurons. In situ hybridization was performed on human embryonic brain sections. Genetically engineered mice were generated with a series of C-terminal OTX2 variants. Results: Two chromosomal deletions and six haploinsufficient mutations were identified in individuals with eye abnormalities; an affected relative of one patient harboured the same mutation without an ocular phenotype. OTX2 truncations led to significant transactivation reduction. A missense variant was identified in another patient without eye abnormalities; however, studies revealed it was most likely not causative. In the mouse, truncations proximal to aa219 caused anophthalmia, while distal truncations and the missense variant were tolerated. During human embryogenesis, OTX2 was expressed in the posterior pituitary, retina, ear, thalamus, choroid plexus, and partially in the hypothalamus, but not in the anterior pituitary. Conclusions: OTX2 mutations are rarely associated with hypopituitarism in isolation without eye abnormalities, and may be variably penetrant, even within the same pedigree. Our data suggest that the endocrine phenotypes in patients with OTX2 mutations are of hypothalamic origin.